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ABSTRACT 

Using the unprecedented spectral resolution of the reflection grating spec- 
trometer (RGS) onboard XMM-Newton, we reveal dynamics of X-ray-emitting 
ejecta in the oxygen-rich supernova remnant Puppis A. The RGS spectrum shows 
prominent K-shell lines, including VII Hea forbidden and resonance, VIII 
Lya, VIII Ly/3, and Ne IX Hea resonance, from an ejecta knot positionally 
coincident with an optical oxygen-rich filament (the so-called "O" filament) in 
the northeast of the remnant. We find that the line centroids are blueshifted 
by 1480±140±60 kins" 1 (the first and second term errors are measurement and 
calibration uncertainties, respectively), which is fully consistent with that of the 
optical Q filament. Line broadening at 654 eV (corresponding to O VIII Lya) is 
obtained to be a < 0.9 eV, indicating an oxygen temperature of < 30keV. Analy- 
sis of XMM-Newton MOS spectra shows an electron temperature of ~0.8keV and 
an ionization timescale of ~2 x 10 10 cm -3 s. We show that the oxygen and electron 
temperatures as well as the ionization timescale can be reconciled if the ejecta 
knot was heated by a collisionless shock whose velocity is ~600-1200 kms -1 and 
was subsequently equilibrated due to Coulomb interactions. The RGS spectrum 
also shows relatively weak K-shell lines of another ejecta feature located near 
the northeastern edge of the remnant, from which we measure redward Doppler 
velocities of 650±70±60 kins -1 . 
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1. Introduction 



The Galactic supernova remnant (SNR) Puppis A is one of the so-called "oxygen-rich" 
SNRs, in which oxygen-rich fast-moving optically emitting ejecta knots (OFMKs) are iden- 
tified. So far, eight other SNRs, Cas A and G292.0+1.8 in our Galaxy and six others in the 
Magellanic Clouds and NGC 4449, are categorized in this class (jVinkll2012l . for a recent sum- 
mary). Of these, Puppis A is particularly i nteresting, because all of OFMKs are located i n 



the eastern and northeastern (NE) portion ( Winkler Sz Kirshner 1985 : Winkler et al. 1988}) 



while the stellar remnant is moving in the opposite southwester n direction ( IPetre et al 



1996; 



Hui fc Becker! 120061 : IWinkler fc Petrel 120071 : Becker et all 120121 ) . This is consistent with re- 
cent models, which predicts that SN ejecta are expelled prefer entially in one direction, while 



the stellar remnant receives a kick in the opposite direction ( iWongwathanarat et al.l 12012 



and references therein). Thus, Puppis A is an extremely important target for the study of 
SN explosion mechanisms as well as SN nucleosynthesis. 

Puppis A is one of the brightest SNRs in the X-ray sk y. The X-ray emission has been 



thoug ht to be dominated by strong cloud-shock interactions (ICharles et al.lll978l ; IPetre et al. 



19821 . references therein). A large area of the remnant was recently mapped by X-ray obser- 
vatories in orbit, i.e., XMM-Newton, Chandra, and Suzaku. They confirmed earlier results 



that the X - ray emission is dom inated by the swept-up interstellar medium (IHwang et al. 



20051 . 120081 ; iKatsuda et al.ll2010l ). On the othe r hand, the data al so revealed signatures of 
SN ejecta in three locations: a Si-rich clum p (IHwang et al.l 120081) . a bright 0-Ne-Mg-(Si- 
S)-rich knot (the ejecta knot called hereafter: Katsuda et al. 20081 ). whose southern portion 
is positionally coincident with an -shaped OFMK (the so-called Q filament named by 



Wink: 



er fc Kirshner 



Katsuda et al. 



1985), and an O-Ne-Mg-rich filament (the ejecta filament called here- 



20101 ). All of them are located in the NE quadrant, further supporting 



after: 

the one-sided ejection of SN debris. 

While picking up ejecta features reveal two dimensional ejecta distributions on the plane 
of the sky, measure ments of Doppler velocities can provide us with l ine-of-sight structures 
(IDewey. et al.ll2010l . for a recent review of kinematics of X-ray SNRs). IKatsuda et al.l (120081 ) 
found that the ejecta knot described above shows blueshifted K-shell lines. The Doppler ve- 
locities were measured to be 3400^gog° km s _1 at the northern knot and 1700^ggo km s^ 1 at the 
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south ern knot which agrees with the optical measurement of ~1500 km s _1 ( Winkler &: Kirshner 
19851 ) . The possible different Doppler velocities between the north and the south led us to 
suggest that the ejecta knot actually consists of two different knots which appear to be 
located close with each other. However, the poor quality of the data prevented us from 
making strong arguments; t he X-ray data were taken by the XMM -Newton European Pho- 
ton Imaging Camera (EPIC: iTurner et al.ll200ll ; IStriider et al.ll200ll ) with moderate spectral 
resolution (E/AE ~ 20), and also the ejecta knot was detected at the edge of the field of 
view (FOV) where calibration uncertainties would be large. 

Here, we report on precise Doppler measurements for the ejecta knot and the ejecta fil- 
ament in Puppis A, based on the unprecedented spectral res olution X-ray spectra ( E/AE ~ 
150) obtained by the Reflection Grating Spectrometer (RGS: lden Herder et al.ll200ll ) onboard 
XMM-Newton. In addition, we give a restrictive upper limit on an oxygen temperature of 
the ejecta knot, for the first time. These new measurements allow us to discuss the nature 
of the ejecta knot in detail. 



2. Observation and Data Reduction 



We performed an XMM-Newton observation of Puppis A on 2012 October 20 (Obs.ID 
0690700201), aiming at the fast-moving ejecta knot coincident with the optical Q filament. 
Our primary instrument in this observation is the RGS that can produce high-resolution 
X-ray spectra, while we also make use of the EPIC data to support the RGS results as well 
as to obtain additional information. The dispersion direction of the RGS is 96°. 4 measured 
east to north, as we show in Figured] left — an X-ray image of Puppis A generated by existing 
XMM-Newton and Chandra data. Fortunately, no background flare is seen in our light curve, 
allowing us to use full exposure times of 20.8ks, 12. Iks, and 11.7ks for the RGS1/2, MOS1, 
and MOS2, respectively. All the raw data were processed using version 11.0.0 of the XMM 
Science Analysis Software and the calibration data files available in 2012 May. 



3. RGS Spectra 

The slitless RGS spectrometer is generally not useful for extended sources like Puppis A, 
because off-axis emission along the dispersion direction is detected at wavelength positions 
shifted with respect to the on-axis source, which causes degradation of energy resolution. 
However, if the angular size of the target is small enough (< a few arc minutes) and is 
brighter than its surroundings, it is possible to obtain high-resolution spectra for them. Our 
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on-axis source, the ejecta knot, is such a lucky case, for which we are able to obtain an order- 
of-magnitude better resolution spectra than nondispersive CCDs. There is another ejecta 
feature, i.e., the ejecta filament described in the introduction, within the RGS FOV. Since 
the angular size of the filament along the dispersion direction is as thin as 1', we can also 
find its signatures in the RGS spectrum. The two ejecta features of interest are indicated on 
a Chandra image shown in Figure [1] right. 



The analysis method basically follows our previous work (IKatsuda et al.ll2012l ). Briefly, 
we first extract four RGS spectra from 0'. 8- width sectors along the cross-dispersion direction. 
We then smooth the RGS response designed for on-axis point sources, based on emission 
profiles along the dispersion direction by using the rgsrmf smooth s oftware. Fina lly, we 



fit the RGS spectrum with emission models in the XSPEC package (lArnaudl Il996l ). The 
RGS background (BG) is generated from a blank-sky observation (Lockman hole: Obs.ID 
0147511601), as we did in our previous work. The blank-sky BG emission is less than 1% of 
the source emission in the energy band analyzed. 

Here, the RGS response in each region is smoothed with three distinct emission profiles 
responsible for (1) the ejecta knot, (2) the ejecta filament, and (3) their surrounding diffuse 
emission, so that we can derive physical parameters for individual features. Emission profiles 
are generated with either the Chandra or XMM-Newton image, for which we take account 
of vignetting effects of XMM-Newton 's X-ray telescope and the asymmetry of the effective 



area along the dispersion axis described in the literature ( IBroersen et al.ll2013f ) . The emission 
profiles for the ejecta features are extracted from Chandra images masked for them, for which 
we subtract underlying diffuse emission by interpolating the surrounding diffuse emission. 
We have checked that, even if we change the intensity of underlying diffuse emission by ±20% 
which is a small-scale inhomogeneity around the ejecta knot, the fit results presented below 
are not affected significantly. Figure H] shows the X-ray emission profiles in 0.6-0.7 keV for 
the four regions. The three profiles in each panel are responsible for (1) the ejecta knot in 
blue, (2) the ejecta filament in red, and (3) the surroundings in black. 

Figures [3] and @] shows BG-subtracted RGS1 and RGS2 spectra for the four regions. 
While most of the emission is due to overall diffuse emission which acts as effective BG 
emission in our analysis, several peaks can be identified as line emission from our interesting 
two ejecta features. Some of them are indicated in the plot as blue and red labels for the 
ejecta knot and the ejecta filament, respectively. Note that there is large wavelength shifts 
in line centroids between the knot (blue) and the filament (red) by about lA. This is mainly 
due to the off-axis position (about 8') of the filament; an offset of 1' causes a wavelength 
shift of 0.138 A. We use data in the energy range of 0.5-1 keV (24-13 A) where line emission 
from the ejecta features are clearly detected. Dividing the energy band into four sections, 
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0.5-0.6 keV, 0.6-0.7keV, 0.7-0.85 keV, and 0.85-0.97keV, we generate RGS responses using 
the same energy-band images. We note that the available energy ranges for the RGS1 and 
the RGS2 are limited to 0.5-0.8 (24-15.5 A) and 0.65-1 keV (19-13 A), respectively, since 
chip #7 in the RGS1 and chip $4 in the RGS2 are no longer operating. 

The RGS spectra consist of diffuse emission and emission from the ejecta knot and 
ejecta filament. We model the diffuse emission an d the ejecta in terms of a nonequilib- 
rium ionization shock model (the vpshock model: iBorkowski et al.l 120011 ) and individual 
Gaussian components, respectively. We use nine gaussian components for the ejecta knot 
and six for the ejecta filament. We impose the Gaussians in the knot component to be 
O VII Hea forbidden at @0.561keV (22.1 A), O VII Hea intercombination at @0.568keV 
(21.8 A), O VII Hea resonance at @0.574keV (21.6 A), O VIII Lyct @0.654keV (19.0 A), 
O VII He/3 @0.666keV (18.6 A), O VIII Ly/3 @0.775keV (16.0 A), Ne IX Hea forbidden 
@0.905keV (13.7 A), Ne IX Hea intercombination @0.915keV (13.6A), and Ne IX Hea res- 
onance @0.922keV (13.5 A). The six Gaussians for the filament component are the same O 
K-shell lines listed above. For the diffuse component, the degradation of RGS spectra is 
sever owing to the large spatial extent; the effective spectral resolution is much worse than 
the nondispersive CCDs. Therefor, it is difficult to constrain numerous parameters, hence 
in our fitting procedure, we only allow the electron temperature, kT e , and normalizations 
to vary freely, fixing the intervening hydrogen column density, Nn, to be 3xl0 21 cm~ 2 and 
the range of the ioni zation timescale, n fi t, from zero up to 3 xlO n cm -3 s, based on recent 



X-ray measurements ( iHwang et al.ll2005l . 12008c iKatsuda et al.ll2010f ). For the Gaussian com- 
ponents, we allo w all normalization s to vary freely. Line centroids are fixed to theoretically 
expected values (ISmith et al.ll200l[ ). but we allow a redshift parameter to vary freely, keep- 
ing the same value for all Gaussians in either the knot or the filament component. Line 
broadening of the knot component is a free parameter taken to be proportional to line center 
energies as is expected for a collisionless shock heating. That of the filament component 
is set to zero, since line emission from the filament is not very clear in the RGS spectrum. 
For the knot component in regions A, B, and C, we measure redshift and line broadening 
separately for the RGS1 and the RGS2, so that we can evaluate systematic uncertainties in 
our measurements and also see if energy dependence is present; the RGS1 covers 0.5-0.8 keV 
including O K-shell lines while the RGS2 covers 0.65-1 keV including O and Ne K-shell lines, 
respectively. Before BG subtraction, each spectrum is grouped into bins with at least ~25 
counts, which allows us to perform a \ 2 test. 

With this fitting strategy, we obtain fairly good fits for all the spectra as shown in 
Figure EJ where individual knot and filament components are shown in blue and red, respec- 
tively, as well as the total model in green. The fit results are summarized in Table [TJ Overall, 
the redshifts and line broadening derived from the RGS1 and the RGS2 are consistent with 
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each other, suggesting no significant energy dependence for these parameters. We also fitted 
the RGS2 spectra, by using a specific energy band for Ne IX Hea lines, finding consistent 
results as in Table [TJ Figure @] shows close-up RGS spectra for O K-shell lines. In this plot, 
we show the best-fit models shifted at rest-frame positions (i.e., redshift = 0) in green. The 
Doppler shifts of the ejecta knot are clearly seen. 

The blueshift measured at the southernmos t region (region D), 1 560±60 kms -1 , is con- 



sistent with our previous X-ray measurements ( jKatsuda et al.ll2008l ). significantly reducing 
the uncertainty. It also matches with the previous optical measurements for the Q fila- 
ment: 1520±40 kms" 1 or 1570±20 kms -1 from the SIT or CCD spectrum, respectively 



( IWinkler fc Kirshnerlll985l ). While the blueshifts in the four regions are all close with each 
other, that in the northernmost region is slightly slower than the other regions. This con- 
flicts with our previous measurements that suggested a faster velocity in the north of the 
knot than that in the south. Given much improved spectral resolution in the current data, 
we believe that the result presented here is more reliable than the previous one. As for 
the slightly slower velocity, it may indicate that the northern part is responsible for a tail 
structure delayed from the main body of the ejecta knot. We find, for the first time, that 
line emission from ejecta filament is redshifted by ~650 kms -1 . 

Line broadening in the ejecta knot is found to be fairly narrow. It ranges from to 0.9 eV 
at 0.654 keV in the four regions, which we take as our measurement unce rtainty. Based on the 



relati on between temperatures and line broadening, a = E ^J kT / mc 2 (IRybicki fc Lightman 



19791 ). we obtain an upper limit of an oxygen temperature to be 30keV. 



Line intensity ratios measured are consistent with those expected by a non-equilibrium 



ioniz ation plasma model. We have indeed confirmed that the vnei model (IBorkowski et al. 
20011 ) can also well reproduce the RGS spectra (values of x 2 /d.o.f. range from 1.4 to 1.6) 



with plausible electron temperatures, ~0.4-lkeV, an ionization timescale of 2xl0 10 cm 3 s, 
and both Doppler shifts and line broadening consistent with those in Table [TJ 



4. MOS Spectra 

Since the RGS spectra are heavily contaminated by the surrounding diffuse emission 
in our case, electron temperatures and ionization timescales for our interesting features are 
better constrained by the MOS rather than the RGS. We thus analyze the MOS spectra for 
the ejecta knot. As shown in a close-up Chandra image in Figure [5j the knot appears to 
consist of three parts: north, west, and south. We extract spectra from these features. Local 
BGs are taken from the adjacent regions. We sum up data taken by the MOS1 and the MOS2 
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to improve photon statistics. The local-BG subtr acted MOS1+2 spectra together with the 
best-fit vnei model (the augmented version 2.0: iBorkowski et al.ll200ll ). are presented in 
Figure El 

It is difficult to measure absolute (relative to hydrogen) metal abundances for metal- 
rich plasmas, because continuum emission comes from the metals themselves and no longer 
reflects the amount of hydrogen, and also because in our case the continuum emission from 
the ejecta knot is difficult to estimate adequately, given that a considerable fraction of the X- 
ray emission is due to local BG that cannot be subtracted perfectly. Therefore, in our spectral 
modeling, we examine two abundance patterns: I n case A, we assume the Fe/H abundance 
to be the solar value ( Anders fe Grevesse I Il989l ) as line emission from Fe is not evident 
in our X-ray spectra, while in case B, we assume the O/H abundance to be 2000 times the 
solar value which is derived from optical spectroscopy of the f2 filament (jwinkler fe Kirshner 
19851 ). We find that the vnei model represents the data fairly well for both of the two cases. 
The best-fit models presented in Figure [6] are those for case A which gives slightly better fits 
than case B. Table [2] summarizes the results and notes details of our fitting. 

The blueshifts from the MOS spectra are in reasonable agreements with the RGS mea- 
surements, considering the relatively large calibration uncertainty (±5eV0) on the energy 
scale of the MOS. We also confirm the elevated abundances of O, Ne, and Mg compared with 
Fe. On the other hand, the electron temperatures a nd ionization timesca les are outside sta- 
tistical uncertainties of our previous measurements (IKatsuda et al.ll2008l ). The discrepancy 
would be partly due to the different data quality; the data were taken at the edge (previous) 
or center (this time) of the FOV, and partly due to the different fitting strategy; we here use 
(1) the most recent version of the vnei code, (2) different spectr al extraction regions , and 
(3) different energy band, and also we fix Nn to be 3xl0 21 cm -2 (IKatsuda et al.l 120101 ). but 
allow abundances of C and N to vary freely. Not only these analysis improvements but also 
the fact that blueshifts derived here are closer to the RGS measurements assure reliability 
of the current results rather than the previous ones. 



Assuming the plasma depths to be 1' (or ~0.6 pc at a distance of 2.2 kpc: iReynoso et al. 



20031 ). which is roughly the size of the knot, we calculate densities and masses in each 
region. Although it is often assumed that a density ratio between electrons and protons to 
be unity or so, this assumption is incorrect for super-metal rich plasmas like the ejecta knot 
in Puppis A. Therefore, we here calcula te the ratio, taking a ccount of abundances as well as 
ion fractions based on the SPEX code (IKaastra. et al.lll996l ). Then, the emission measures 
and abundances give densities and masses for each species as summarized in Table [31 The 



1 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.ps.gz 
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total mass in each region, ~0.01 M , is consistent with our previo us estimates ( Katsuda et al. 



20081 ) . and is comparable with typical OFMKs in the remnant (IWinkler et al. 



1988). 



Discussion 



The XMM-Newton RGS allowed successful measurements of accurate Doppler velocities 
of two X-ray-emitting ejecta features (i.e., the ejecta knot and the ejecta filament) in the 
NE quadrant of Puppis A. In our spectral analysis, we divided each feature into four cross- 
dispersion regions. The error-weighted mean Doppler velocities for the four regions are 
obtained to be 1480±140±60 kms -1 blueward for the ejecta knot and 650±70±60 kms -1 
redward for the ejecta filament, where the first and second term errors, respectively, represent 
the standard deviation for the four regions and a wavelengths accuracy of the RGSo The 
Doppler velocity o f the ejecta knot is fully c onsistent with previous optical measurements 
for the Vt filament ( Winkler fc Kirshnerlll985l ). ensuring robustness of our measurements. 



The RGS spectra also enabled, for the first time, to measure broadening of O (and 
Ne) K-shell lines of the ejecta knot. The broadening is found to be modest, a < 0.9 eV 
at O VIII Lyct, indicating an upper limit of an oxygen temperature of 30keV. As far as 
we know, this is the second result that we were able to measure an oxygen temperature of 
fast-moving ejecta knots in X-rays. Comparing our results with the first successful object, a 
comparably fast-moving ejecta knot in the northwestern rim of SN 1006, we notice that the 
oxygen temperature in Puppis A's knot is at least an order-of-magnitu de lower than that in 
SN 1006's knot (kT Q ~300keV: Ivink et alJliooit iBroersen et aDl2013h . 



The post-shock temperature, Tj, for species i with mass rr ij right behin d a collisionless 



shock in the absence of the magnetic pressure is given as (e.g., ISpitzerlll965l ): 



kT, 



16 



(1) 



for shock velocity, t> s h- Thus, the particles will have mass-proportional temperatures in the 
limit of collisionless plasma. In fact, the high oxygen temperature together with a relatively 
low electron temperature of ~1.5keV in SN 1006 has been reasonably interpreted as evidence 
for temperature nonequilibration between ions and electrons after collisionless shock heating. 
As for Puppis A's knot, the oxygen temperature is expected to be ~130keV, if we assume 
that the shock velocity is similar to the gas velocity: ~ 2000 kms _1 derived by the Doppler 
velocity of ~ 1500 kms~ and the optical proper motion of ~1250 kms -1 at a distance 



2 http: / / xmm2.esac.esa.int / docs / documents / CAL-TN-0030.pdf 
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of 2.2 kpc flWinkler et all 11988 : iGarber et al.l I201CH : iRevnoso et all 120031 ). This prediction 
conflicts with our measurement, kTo <30keV. 

To understand the apparent temperature inconsistency between the measurement and 
the simple prediction above, we are required to inspect heating processes at the collisionless 
shocks. In fact, it is well known tha t temperatures measured in the im mediate post-shock 
regions do not often obey Equation [Tl (jGhavamian et al.ll2007l : IVinkll2012L for a recent review). 
The deviation may be due to plasma processes that ra ise an electron temperature at (in front 
of) the shock (jOhira et al.ll2008l ; iRakowski et al.ll2008l ) and/or significant energy leakage into 



cosmic rays that l owers post-shock temperatures of all particle species ( [Hughes et al.l 12000 



Helder et al.l 120091 ). In any case, oxygen can have a lower temperature than that predicted 
by Equation [TJ 

In this context, it is important to clarify whether the Tq/T s ratio observed in Puppis A's 
knot can be explained by temperature equilibration due to Coulomb interactions behind a 
collisionless shock or requires other scenarios. To this end, assuming pure collisionless shock 
heating, i.e., initial temperatures are given by E quation HQ we s olve the following temperature 
equilibration due to Coulomb interactions (e.g., SpitzerT 1965 ): 



dt 



C eq(i,j) 



C eq(i,j) 



5.87 



AiAj 



. 3 
TA* 



T; 



'-J 



(2) 



for species i and j with atomic numbers A4 and Aj, charges and Zj, number density rij, 
and the Coulomb logarithm ln(A) ~ 40. We coupled the differential equations for electrons, 
protons, He, C, O, Ne, Mg, Si, and Fe, using the abundances for either case A (He/H=l, 
C/H=0/H=5, Ne/H=8, Mg/H=7, Si/H=3, and Fe/H=l times the solar values, see Ta- 
ble [2D or case B (He/H=l, C/H=O/H=2000, Ne/H=3000, Mg/H=2500, Si/H=1000, and 
Fe/H=300 times the solar values). After examining several shock velocities, we find that the 
electron temperatures and ionization timescales in the ejecta knot can be best explained at 
v s h ~ 1200 kms -1 or v s h ~ 600 kms -1 for case A or case B, respectively. The temperature 
histories for electrons, protons, and oxygen as a function of n e t are shown in Figure [7J From 
the plot, we also find that the oxygen temperature measured, T G < 30 keV, is consistent with 
the expectation in both cases. Therefore, we do not need to invoke some exotic scenarios to 
explain the low To/T e ratio. 

Still, the temperatures measured may be also reproduced for the case that initial post- 
shock temperatures are not proportional to particle masses, as is often observed. If it is 
the case, the shock velocities would be different, causing uncertainties on the shock veloc- 
ity. Recent observational studies agree in that the initial electron-to-proton temperature 
ratio, T e /Tp = (3, is close to unity (full equilibration) for slow shocks (v s h <400 kms -1 ), 
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but gradually goes down for faster shocks (jGhavamian et al.l 120071 ; iHend |2010| . and refer- 
ences therein). We take th e value of f3 to be 0.2 which is typical around t> s h = 1000 kins -1 
(Ivan Adelsberg et al.ll2008l ). We find that temperature histories expected at t> sh ~ 1000 km s -1 
(case A) or v s h ~ 600 kins -1 (case B), indicated as dotted lines in Figure [3, best match our 
measurements. In this context, we argue that the shock velocity does not depend on the 
degree of initial temperature equilibration compared with the ambiguity of abundances. 

It is interesting to note that the inferred shock velocity, 600-1200 kms -1 , is less than 
the shocked-gas velocity, ~ 2000 kms -1 . This suggests that the ejecta knot was heated 
by a reverse shock rather than a forward shock, since the forward shock velocity should be 
~2500 km s -1 according to the Rankine-Hugoniot relation: 



7 + 1 



(3) 



where Vf s and v sg are velocities of the forward shock and the shocked gas, respectively, and 7 is 
the specific heat ratio taken to be 5/3 for a non-relativistic ideal gas. The Rankine-Hugoniot 
relation over the reverse shock can be written down as: 

7 + 1 / N 



^free 



7-1 



(4) 



where ff ree and t> rs are velocities of the freely-expanding ejecta and the reverse shock in the 
observer's frame. Substituting f free — t> rs (=f s h) = 600-1200 kms -1 , v sg = 2000 kms -1 , and 
7 = 5/3, we compute the value of t>f ree to be 2900-2450 kms -1 , depending on the shock 
velocity of 600-1200 k ms" 1 . This velocity agrees with a typical eject a velocity in O-rich 
layer for Type lib SN (IHouck fc Franssonl Il996l ; llwamoto et al.l 119971 ) which is likely the 



origin of Puppis A (I Chevalier 



20051 ). Also, the distance, 0.3-0.6 pc, that the shock of this 



velocity travels for ~500 years (which is a product of the ionization timescale, 3 x 10 cm s, 
divided by the electron density, 2 cm -3 , as in Tables [2] and [3]) is comparable with the size of 
the knot (~0.6pc), showing self-consistency of our interpretation. 

We comment on the remarkable difference in Tq/T s between Puppis A's knot and 
SN 1006's knot. The proper mo tion of SN 1006's kn ot is measured from recent X-ray 
observations to be ~3000 kms -1 (IKatsuda et al.ll2013l ). Provided that it is shocked by a 
forward shock, then the shock velocity would be the proper motion itself. On the other 
hand, if it is shocked by a reverse shock, the proper motion should be the velocity of the 
shock ed gas. Also, freely-expanding ejecta velocities is likely ~10000 kms -1 (llwamoto et al. 



1999|), since there is a general consensus that SN 1006 is a Type la SNR. This is indeed 
comparable to a freel y-expanding Si-rich e jecta velocity of ~7000 kms -1 measured by the 
UV absorption edge ([Hamilton et al.l 120071 ) as well as a simple estimate from the radius at 
the knot (~9pc) and the SNR age (~1000yr). Then, Equation H] gives the reverse shock 
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velocity in the ejecta-rest frame to be ~9300 kms~ (or ~5300 kms -1 , if we take t>f ree to be 
7000 kms -1 ). Therefore, the shock velocities estimated for SN 1006's knot are at least a 
few times faster than that for Puppis A' knot. In addition, the ionization timescale is an 
order-of-magnitude lower in SN 1006's knot than in Puppis A' knot. These differences can 
qualitatively explain the different oxygen temperatures. 

Future high-resolution X-ray spe ctroscopy with the nondispersive soft X-ray spectrom- 
eter onboard Astro-H (FWHM~5eV: iTakahashi et al.lboid iMitsuda et aDboioh will reveal 
ion and electron temperatures for many SNRs. Hopefully, we will be able to measure ion 
temperatures in different ionization states for every species. Such information will signifi- 
cantly improve our understanding on temperature equilibration for fast collisionless shocks. 
One particularly interes ting; target in Puppis A would be the Si-rich clump in the NE rem- 
nant (jHwang et al.ll2008f ). for which the RGS cannot produce high-resolution spectra due to 
its large angular size (3' x 5'). Fortunately, this target is scheduled to be observed in 2013 
January with a soft X-ray spectrometer onbo ard the Micro-X sounding rocket experiment^] 
(FWHM~2eV: iFigueroa-Feliciano et al.ll2012l ). This experiment will provide exciting results 
shortly. 



6. Summary 

We have presented results from our new XMM-Newton RGS and MOS observations 
of two ejecta features in the Galactic SNR Puppis A: one is a n ejecta knot positionally 
coincident with the optical Q filament (I Winkler fc Kirshnerlll985l ). and the other is an ejecta 
filament located near the NE edge of the remnant. The results obtained are summarized 
below: 



• The Doppler velocity of the ejecta knot is measured to be 1480±140±60 kms -1 blue- 
ward, which is fully consist ent with previous optical measurements for the Q filament 
(IWinkler &: Kirshnerlll985l ). We find that the Doppler velocities are constant within 
the knot, except for a northernmost piece where a slightly slower velocity is indicated. 

• An upper limit of line broadening for O K-shell lines in the ejecta knot is obtained to be 
~0.9eV, leading to an oxygen temperature of < 30keV. This upper limit, combined 
with the electron temperature and the ionization timescale derived from the MOS 
spectra, suggests that the ejecta knot was heated by a reverse shock whose velocity is 



3 http: / /space. mit.edu/micro-x/index. html 
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~600 1200 kms , without requiring rapid temperature equilibration due to plasma 
instabilities nor temperature dropping due to energy leakage into cosmic rays. 

• The oxygen temperature of the ejecta knot is an order-of-m agnitude lower than that ob- 
taine d in a comparably fast-moving ejecta knot in SN 1006 fjVink et al.ll2003l ; iBroersen et al 



20131 ). We argue that the difference would be due to both different shock velocities 



and the ionization timescale. 



• The Doppler velocity for the ejecta filament is also measured, for the first time, to be 
650±70±60 kms -1 redward. 



We thank the referee, Knox Long, for a number of comments that help clarify the paper. 
S.K. is supported by the Special Postdoctoral Researchers Program in RIKEN. This work 
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Table 1. Spectral- fit paramters derived from the RGS 



Com pone nt 


Parameter 




Regions 










A 


B 


C 


D 


Knot 


Rcdshift (10~ 3 ) 


_4 5+ u - a _4 i+ (, -l 


-5.0±0.2, -4.6±0.2 


os -o.i' ~ J - J -0.1 


-5.2±0.2 




Line broadening (<x in cV <§ 


10.654 keV) 0.0 (< 0.4), 0.0 (< 0.7) 


"•3lo.2> 0-9±0-3 


0.0 (< 0.4), 0.0 (< 0.4) 


0.8±0.3 




Intensities: O Hcck (f) 


411±101 


812±123 


882±126 


657±104 




O Hco (i) 


6 (< 87) 


(< 49) 


113±104 


(< 74) 




O Hect (r) 


702±103 


1623±133 


1519±127 


1085±116 




O Hc/3 


157±30 


342±38 


428±40 


134±33 




O Lya 


956±54 


1898±64 


863±52 


220±37 




O Ly/3 


123±20 


247±24 


107±21 


40±16 




No Hca (f) 


26±25 


115±29 


139±28 


48±22 




No Hca (i) 


(< 6) 


25 (< 54) 


42±28 


9 (< 30) 




Nc Heo (r) 


149±25 


216±29 


191±28 


47±22 


Filament 


Rcdshift (10~ 3 ) 


2 ,+0.ii 
"-0.5 


1 9+ U4 


2 2+ U1 


2 4+°- 4 




Intensities: O Heoj (f) 


285±112 


489±118 


374±123 


207±179 




O Hect (i) 


(< 48) 


77 (< 185) 


70 (< 167) 


129±104 




O Hco (r) 


459+101 


612±109 


1069±119 


619±117 




O He/3 


26 (< 63) 


28 (< 69) 


213±38 


366±43 




O Lya 


237±41 


167±40 


210±38 


223±40 




O Lyp 


113±35 


56±34 


(< 34) 


« 7) 


Diffuse 


kT (kcV) 


0.46±0.01 


0.47±0.01 


0.44±0.01 


0.45±0.01 






1937.8 / 1260 


2073.0 / 1280 


1841.6 / 1260 


1875.9 / 1257 



Note. — Wh values are fixed to 3xl0 21 cm -2 . Line intensities arc in units of 10 -5 photons cm -2 s~ 1 , respectively. Two values of the rcdshifts 
and line broadening in Regions A, B, and C are derived by the RGS1 (left) and the RGS2 (right). The errors represent 90% confidence levels on an 
interesting single parameter. 



Table 2. Spectral-fit parameters derived from the MOS 



Parameter North West South 





Case A 


Case B 


Case A 


Case B 


Case A 


Case B 


kT c (kcV) 
log(n c t /cm ^ s) 


o.93+. ;i b 2 
io.43± u ; o 6 


1.03+0.12 
10 40+ 07 


80 + ,, UU 

10.24+g;S? 


S3+ U lj(i 
0.83_ 12 

10 23+ 07 


,. ,..,+o.ia 

10 21+ 09 


63+ U1(i 
10 22+ 09 


(C/H)/(C/H) 


11 9+ 80 
iJ -' J -7.0 


KMMM+2400 




1800+1800 


3 1+ 36 


1500+1300 


(N/H)/(N/H) 


0.0 (< 1.2) 


(< 300) 


1.8 (< 3.8) 


500 (< 1100) 


1 5+ 1 - 4 


700+500 


(O/H)/(O/H) 


, ,+ 1.1 
°-^-0.8 


2000 a 


6 9+ 11 
°- y _0.8 


2000 a 


4 5+ 1 - 1 


2000 a 


(Ne/H)/(Nc/H) 


8.5ti;t 


2400+ 100 


10.3+2- 4 


3000+200 


7 g + 2.2 
'•°-1.7 


3400+300 


(Mg/H)/(Mg/H) 


5 6+ 1 

O.D_ g 


2100+200 


79 +0.7 
l,s -1.2 


2300+200 


7.2+1-9 


3200+500 


(Si/H)/(Si/H) 


3.2+0.7 


1200t?gg 


2.8_ x 2 


900+300 


4 0+ 1 - 9 


1900+ 600 

1MUU_ 500 


(Fc/H)/(Fc/H) 
Redshift (10~ 3 ) 
n n H V/47rd 2 (10 7 c m - 5 ) 


l a 

6 5+ 06 
740 + 130 


380t™ 
-6 6+ 1 - 
1 84+ - 17 


l a 
-6 1+ 03 
240+.™ 


270+ 90 
-6 l+° 3 


l a 

-5 4+ 1 - 6 
°*-0.2 

460+ 170 


390+ 130 
-5 4+ 1 - 
1 04+ - 15 


xVd-o.f. 


187.7 / 126 


200.6 / 126 


141.0 / 103 


148.3 / 103 


180.8 / 108 


192.1 / 108 



Note. — a (Fc/H)/(Fe/H) Q and (O/H)/(O/H) arc fixed to 1 (Case A) and 2000 (Case B), respectively. N H values 
are fixed to 3xl0 21 cm~ 2 . The errors represent 90% confidence levels. 
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Table 3. Densities (10 4 cm 3 ) and Masses (10 4 M Q ) in the ejecta knot 



Parameter North West South 





Case A 


Case B 


Case A 


Case B 


Case A 


Case B 


n G 


44000 


11000 


33000 


7800 


37000 


7200 


"H 


35000 


360 


26000 


370 


29000 


340 


no 


150 


T70 


53 


240 


33 


100 


no 


150 


610 


150 


620 


110 


570 


"•No 


28 


110 


33 


130 


27 


140 


»M| 


7.5 


28 


7.7 


31 


8.1 


41 


"Si 


4.0 


15 


2.6 


11 


4.2 


23 


"Fc 


1.6 


6.3 


1.2 


4.7 


1.4 


6.2 


M H 


82 


0.8 


36 


0.5 


62 


0.7 


M C 


4.3 


22 


14 


4.0 


0.8 


4.7 


Mo 


5.8 


23 


3.3 


14 


3.8 


19 


*f No 


1.3 


5.0 


0.9 


3.7 


1.1 


5.9 


M Mg 


0.4 


1.6 


0.3 


1.0 


0.4 


2.1 


M Si 


0.3 


1.0 


0.1 


0.4 


0.2 


1.3 


M Fc 


0.2 


0.8 


0.1 


0.4 


0.2 


0.7 


Note. - 


- Typical errors 


which arc 


mainly due to the at 


isumption c 


■f the plasma 


depth arc 


about a factor 


of 2. 












Fig. 1. — Left: The RGS spectral extraction regions overlaid on the merged XMM-Newton 
and Chandra image in the energy band 0.5-5 keV. The arrow indicates zero points of projec- 
tion profiles shown in Figure [2J Right: Three-color close-up Chandra image of the white box 
region in the left panel. Red, green, and blue correspond to 0.5-0.7keV, 0.7-1.2 keV, and 
1.2-5.0 keV, respectively. Two ejecta features focused in this paper are indicated as arrows. 
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Fig. 2. — X-ray emission profiles in 0.6-0.7keV along the RGS dispersion axis for the four 
spectral extraction regions. The three profiles in each panel are responsible for the ejecta 
knot in blue, the ejecta filament in red, and the diffuse surroundings in black. The X- 
ray intensities are normalized at the peak values of the ejecta knot, showing how diffuse 
surrounding emission is significant in the RGS spectra. The zero points in x-axis (x=0) 
correspond to the direction from which emission is detected at the nominal (no red/blue 
shift) wavelength positions on the RGS detectors. 
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Fig. 3. — XMM-Newton spectra of the four regions in Figure [U The RGS1 and RGS2 cover 
15.5-24 A (~0.8-0.51 keV) and 13-19.1 A (~l-0.65keV), respectively. The two (RGS1 and 
RGS2) data sets are simultaneously fitted with a combination model of vpshock (for diffuse 
background emission) + nine Gaussians (for the ejecta knot) + six Gaussians (for the ejecta 
filament). The total best-fit model is shown in green, while individual components are in 
blue and red for the ejecta knot and the ejecta filament, respectively. Lower panels show the 
residuals. 
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Fig. 4. — Close-up RGS spectra for K-shell lines. To clarify the Doppler shifts of the two 
ejecta features, we here show the best-fit models shifted at rest- frame positions, i.e., redshift 
= 0. The lower panels show residuals. Blueshifts for O Lya and O Hea for the ejecta knot 
component are clearly seen. 



Fig. 5. — Spectral extraction regions for the MOS data overlaid on a close-up Chandra image 
shown in Figured] right. The local BG regions are also shown. 
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Fig. 6. — The MOS spectra extracted from the regions shown in Figure [5j The best-fit vnei 
models are shown together in green. For clarity, the spectra from the north and the south 
are multiplied by 5 and 0.2 times, respectively. The lower panels are the residuals for north, 
west, and south from top to bottom, respectively. 
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Fig. 7. — Temperature histories for electrons (black), protons (red), and oxygen (blue) as 
a function of n e t. The models are calculated for two initial electron-to-proton temperature 
ratios, 0, of the mass ratio (solid curve) and 0.2 (dotted curve). The crosses are electron 
temperatures and ionization timescales measured for the three regions in the ejecta knot 
(see, Table |2]). Left and right panels are calculated for case A (several-solar abundances as 
indicated on top of the plot) and case B (super-solar abundances), respectively. 



